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Recently Mayer et al. [Phys. Rev. E 68, 016116 (2003)] proposed a new way to compute numeri- 
cally the fluctuation-dissipation ratios in nonequilibrium critical systems. Using well-known facts of 
nonequilibrium critical dynamics I show that the leading contributions of the quantities they con- 
sider are in fact one-time quantities which are independent of the waiting time. The ratio of these 
one-time quantities determines the slope of the straight lines observed in the fluctuation-dissipation 
plots of Mayer et al. 



In a recent work [1] Mayer, Berthier, Garrahan, and 
Sollich (MBGS in the following) presented a study of 
ageing phenomena taking place in nonequilibrium Ising 
models in one and two space dimensions after a quench 
from infinite temperature to the critical point located 
at T = T c . MBGS present inter alia Monte Carlo 
simulations at the critical point of the two-dimensional 
Ising model. These simulations are aimed at computing 
fluctuation-dissipation ratios X(t,t w ) which are defined 
by 



X(t, t w ) — T c Rk=o(t, t w )/ 



dCk=o{t,t w ) 
dt w 



(1) 



where t is the time elapsed since the quench (called obser- 
vation time) and t w < t is the waiting time. Ck=o(t,t w ) 
and Rk=o(t,t w ) are the long- wave-limits of the Fourier 
transforms of the commonly studied spin-spin-correlation 
function and of the conjugate response function [2-4]. 
The quantities dCk =°( t - t ^) anc j R^—Q{t,t w ) have been in- 
vestigated field-theoretically by Calabrese and Gambassi 
in [5] . From general scaling arguments they are expected 
to scale in the ageing limit t w 3> 1, t — t w 3> 1 as [5] 

dc k = {t,t w )_ = Agc{t _ tw)a rt\° Fgc{tw/t) (2) 



dt„ 



Rk=o{t,tw) — A R (t — t w ) a 



F R (t w /t) 



(3) 



with a + 1 = ^ and 6 = ^ - ^ - a. Here d is the 

z z z 

number of space dimensions, z the dynamical exponent, 
A c the autocorrelation exponent, whereas r\ is the usual 
equilibrium critical exponent. The functions Fgc(v) and 
F R (v) are universal with F dc {0) = F R (0) = 1. Eqs. (2) 
and (3) can also be written in the following form: 



d Ck=p(t, t w ) 
dt w 



= C fdc(t/t w ) 



Rk=o(t, t w ) — t^ f R (t/t w ) 
where the scaling functions fdc( x ) and f R (x) vary as 



(4) 
(5) 



for x > 1, i.e., 1 < t w < t. Here 9' = a + 9 is the 
well-known initial-slip exponent of the magnetization [6] 
which for the critical Ising model takes the value 0.19 in 
two dimensions. This power-law behaviour (6) will be of 
importance in the following. 

In their simulations MBGS do not have direct access 
to Rk=o{t,t w ) and dCk =°^ t ' t ^ _ They instead investigate 
integrated quantities: 



G(t,t w )= I du 9 Ck ° = C k =o{t, t) - C k =o(t, t w ) 



du 



(7) 



and 



Xm(t,t w )=T c / duR k=0 (t,u) 



(8) 



fdC,B,( X ) ~ x ° 



(0) 



as these quantities are easily obtained in numerical sim- 
ulations. Plotting Xm(t,t w ) against G(t,t w ) they ob- 
tain within the accuracy of their numerical data straight 
lines with constant slopes. The value of the slope 
is identified by MGBS with the fluctuation-dissipation 
ratio X(t,t w ), see Eq. (1), yielding the claim that 
X(t,t w ) is independent of the waiting time t w . Note 
that this independence on the waiting time is not sup- 
ported by the field-theoretical results [5] which yield 
a fluctuation-dissipation ratio (1) dependent on t w at 
two loops. In addition, MBGS conclude that their 
ratio Xm(t,t w )/G(t,t w ) gives the limit value X°° = 

lim ( lim X(t,t w )) for all times t. 

It is the purpose of this Comment to discuss the inte- 
grated quantities involved in the MBGS analysis of the 
critical two-dimensional Ising model. I shall show that 
the leading contributions to (7) and (8) do in fact not 
depend on the waiting time (these one-time quantities 
will be called non-ageing in the following). Furthermore, 
I shall demonstrate that the constant slope observed by 
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MBGS is given by the ratio of the waiting time indepen- 
dent quantities. It is therefore not a direct manifestation 
of ageing, i.e. waiting time dependent, behaviour. 

The origin of the leading, waiting time independent 
term is readily understood by looking at the integrals 
(7) and (8). Inserting the scaling forms (4) and (5) 
one obtains: G{t,t w ) ~ t^ +1 f G (t/t w ) and Xm{t,t w ) ~ 
1 /x^Aw)- However in doing so we did not pay atten- 
tion to the conditions of validity of (4) and (5). Indeed, 
close to the upper integration limit these scaling forms 
cannot be used, since there the condition t 3> t w 3> 1 is 
not fulfilled. One might therefore argue that a time scale 
t* exists such that only for t w < t* the forms (4) and 
(5) hold. As shown in the following, the time integrals 
in (7) and (8) indeed yield a contribution with a scaling 
behaviour which differs from that of an ageing quantity. 

The correlation C).—o(t, t w ) of the magnetization is 
given by [7] 
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C k = (t,t w ) = N M-^«j(<) 




(9) 



where Sj(t) is the value of the spin located at the lat- 
tice site i at time t. N is the total number of lattice 
sites, whereas (• • •) indicates an average over the ther- 
mal noise. In G(t,t w ), analyzed by MBGS, this age- 
ing quantity is subtracted from the one-time quantity 

C k=0 {t,t) = N ((j? V see ( 7 )- This latter 

quantity is usually denoted by M^ 2 \t) in the literature 
and has been extensively studied in the context of short- 
time critical dynamics, see, e.g., [8-10]. Standard scaling 
arguments [9] show that immediately after the quench 
Ck=o(t,t) grows as £( d ~ 2 ^ '< >> lz where d is the number 
of space dimensions and (3 and v are the usual equilib- 
rium critical exponents. For the two-dimensional Ising 
model we have (d - 20/v)/z w 0.81. Therefore G(t,t w ) 
is composed of a non-ageing part (i.e. a part which does 
not depend on the waiting time) and of an ageing part 
where the first one grows much faster in time than the 
second one. Indeed, it follows from the dynamical scal- 
ing behaviour [5] that for later times Ck=o(t,t w ) ~ t e 
with 6' — (d — X c )/z. Rigorous arguments [11] yield the 
inequality A c > d/2, thus that Ck=o(t,t w ) never grows 
faster than t d / 2z . This is illustrated in Fig. 1 where I 
plot G(t, t w ) as a function of t — t w for two of the waiting 
times (t w = 46, 193) considered by MBGS and compare 
them to Ck=o{t,t) (grey line). These data have been 
obtained in standard simulations of the two-dimensional 
Glauber-Ising model with heat-bath dynamics. The two 
dot-dashed lines indicate the two different power laws 
involved: the leading contribution ~ i 81 and the sub- 
leading, ageing contribution ~ t 019 . Clearly G(t,t w ) in- 
creases much faster then expected for an ageing quantity 
and rapidly displays a behaviour similar to its leading 



contribution Ck=o{t, t). As shown in the inset Ck=a(t, t w ) 
itself indeed increases as (t/t w )°- 19 (dot-dashed lines), in 
complete agreement with the general scaling arguments 
given above. This power-law behaviour is already en- 
countered for observation times slightly larger than the 
waiting time. 




FIG. 1. G(t,t w ) vs t — t w for two different waiting times 
t-w = 46 and 193. The grey line is the leading waiting time 
independent contribution Ck=o(t, t) which grows in time with 
a power-law with an exponent (d — 2j3/v)/z ~ 0.81. Systems 
with 300 x 300 spins have been simulated, and all the data 
shown in the present work have been averaged over 70000 dif- 
ferent runs. The inset shows that the ageing part Ck=o(t,t w ) 



of G(t,t w ) indeed grows as t 
guments. 



as predicted by scaling ar- 



To compute the susceptibility of the magnetization 
a small homogeneous constant field of strength h is 
switched on after the waiting time t w [1] . The integrated 
susceptibility is then given by [7] 



N 



Xm(t,t w ) = -^ri^Si) = 



(10) 



The magnetization m is measured for times t > t w and 
depends both on t and t w . Starting from an uncorrelated 
initial state, the dynamical correlation length increases 
with time, £(t) ~ t 1 ^, up to t w . The homogeneous ex- 
ternal field, which is applied for times t >t w , drives the 
system away from the critical point towards a new equi- 
librium point located at T = T c and m — mfi na i > 0. 
This new equilibrium point is reached at finite times, in- 
dependently of the waiting time. This is also the case 
when the system is already in equilibrium at the critical 
point and an homogeneous external field is then switched 
on. One therefore expects that the extension of the cor- 
related regions (and therefore the value of t w ) is only of 
importance for a short period after the application of the 
field, but that at later times the system looses the mem- 
ory of the value of £(t w ), thus that Xm(t, tw) approaches 
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XFc{t) ■= Xm(t,0) for all waiting times. Here, XFcif) 
is the field cooling susceptibility. This is illustrated in 
Fig. 2 where the field strength h = 0.0004 is the same 
as that used by MBGS [12]. One indeed observes that 
the curves with t w ^ rapidly approach the curve for 
t w = 0. Note that the plateau reached at longer times 
is not a finite-size effect, but is due to the new equilib- 
rium point. One further remarks from Fig. 2 that XFc(t) 
exhibits a power-law behaviour. The value of the corre- 
sponding exponent can be obtained from the standard 
dynamical scaling relation for m [9] (with m(t = 0) = 0) 

m(t, r, h) = b- & l v m(b- z t, b 1 ^ t, b d ~ & l v h) (11) 

where r is the reduced temperature. In our case r = as 
the temperature is fixed at T c after the quench. Setting 
b = t x l z one gets 



m(t, h) = r p l» z m(l, t^-P'^l* h) 



and 



T 

XFc(t) = ym(f,/i) 



t (d-2/3/u)/z 



(12) 



(13) 



where the last step h is valid for t < h^ (d ^ f3/l,s>/z . This 
expected power-law behaviour is also shown in Fig. 2. 
It is important to note that XFc(t) increases with the 
same exponent (d — 2(3 jv) / z as the leading contribution 
Ck=o{t, t) of the correlation G(t, t w ), and this already af- 
ter a few time steps. 
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FIG. 2. The susceptibilty of the magnetization x™(~t, t-w) as 
function of the time t — t w elapsed since the switching on of 
the homogeneous constant external field. One observes that 
the curves for waiting times t w 7^ rapidly converge to the 
curve obtained for t w = where the system is quenched to 
Tc in presence of the field. The leading contribution of the 
susceptibility also grows with a power-law with an exponent 
(d — 2(5/v)/z ~ 0.81. The thin dashed line indicates the final 
value of Xm(t,t w ) at the new equilibrium point. 



As the leading non-ageing contributions of both quan- 
tities used by MBGS display the same power-law be- 
haviour, one may wonder whether the straight lines ob- 
served in their fluctuation-dissipation plots, Figs. 8 and 
10 in [1], are not simply due to this non-ageing parts. 
From the scaling arguments one expects that the ra- 
tio Xm(t,t w )/G(t,t w ) is given by Xm{%t w )/G{t,t w ) = 
XFc(t)/C k=0 (t, t) + 0{fP/«'-^l*) with 2(3 /vz - \ c /z = 
—0.62 for the two-dimensional Ising model. The ratio 
XFc(t)/Ck=o(t, t) is expected to take a constant value X 
already after a few time steps. This is indeed the case, 
as shown in Fig. 3. Here I plot XFc(t) as function of 
Ck=o {t, t) and compare the resulting line with those ob- 
tained when plotting Xm(t, t w ) as function of G(t, t w ) for 
t w = 46 and 193, as done by MBGS. The behaviour of 
these quantities at the very first time steps is displayed in 
the inset. It is now obvious why MBGS obtain straight 
lines with a constant slope for any value of the waiting 
time t w : the ageing (i.e. waiting time dependent) parts 
are rapidly suppressed in this kind of plot and the slope 
is then identical to the slope obtained from two quanti- 
ties which do not depend on the waiting time and which 
furthermore have the same time dependence. 
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FIG. 3. Fluctuation-dissipation plot similar to Fig. 8 in [1] 
for waiting times t w = 46 and 193. Only every 200th data 
point is shown. The grey line is the correponding curve ob- 
tained from plotting the leading waiting time independent 
contribution \Fc(t) = Xm(t,0) of the susceptibility as func- 
tion of the leading waiting time independent contribution 
Ck=o(t,t) of G(t,t w ). Obviously, the slopes of the different 
curves are identical after very few time steps and are exclu- 
sively due to the non-ageing parts of the different quantities, 
see inset. 

In conclusion I have shown that the leading terms of 
the integrated quantities used in [1] for the numerical de- 
termination of the fluctuation-dissipation ratios are inde- 
pendent of the waiting time. I have also shown that the 
leading terms of the correlation (9) and the susceptibility 
(10) grow in time with the same power-law. This explains 
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why MBGS observe in their fluctuation-dissipation plots 
straight lines with a constant slope that does neither de- 
pend on the waiting time t w nor on the observation time 
t. 
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